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We investigate a cylindrical plasma system with ionization, by an annular beam of ultraviolet light, taking 
place only in the cylinder’s outer region. In the steady state both the outer and inner regions contain a 
plasma, with that in the inner region being uniform and field-free. At the interface between the two regions 
there is an infinitesimal jump in ion density, the magnitude approaching zero in the quasi-neutral (A^i —> 0) 
limit. The system offers the possibility of producing a uniform stationary plasma in the laboratory, hitherto 
obtained only with thermally produced alkali plasmas. 


I. INTRODUCTION 

In two recent papers we have presented models of a 
plasma split into two regions of equal width, with ion¬ 
ization by ultraviolet light taking place only in the inner 
region and not in the outer— At the interface between 
the two regions a double layer was formed, the ion ve¬ 
locity entering the outer plasma being greater than the 
Bohrn (or ion-acoustic) velocity. In the present paper we 
study the inverse case, i.e. ionization due to photo-violet 
irradiation only in the outer region. A uniform plasma 
fills the inner region, but instead of a double layer at 
the interface between the regions there is an infinitesi¬ 
mal jump in ion density; quasi-neutrality then holds up 
to the formation of a sheath at the wall. In this case the 
ion velocity on entering the wall sheath is about equal to 
the Bohm velocity. The model predicts a uniform station¬ 
ary plasma in the inner region. Such a state of plasma 
is extremely rare; the only case known to us is that of 
a thermally produced alkali plasma-. We can note that 
these two different situations have the common feature 
that the ions and electrons, which constitute the plasma, 
have entered the plasma from outside. The system offers 
the possibility of interesting experimental work, since it 
is difficult to make measurements on thermally produced 
plasmas. 

II. THEORETICAL MODEL 

Figure [1] illustrates the plasma modelled here, which 
four steady-state equations describe. Two are fluid equa¬ 
tions, namely the continuity equation 

V • (rijv) = g (1) 

and the ion momentum equation 

miriiiy ■ V)v + rriiVg = — rijeW (2) 


wall 



FIG. 1. The cylindrical system we analyze. The system com¬ 
prises two regions of plasma, the central region I (0 < r < a) 
and the outer region II (a < r < 2a). Generation by pho¬ 
toionization occurs in region II with intensity gn, and not at 
all in region I. 

while the other two equations are Poisson’s equation 

V 2 y=— (ne-rh) (3) 

£o 

and the Boltzmann relation for electrons 

= B “ exp (i£rj (4) 

where ni and n e are the ion and electron densities respec¬ 
tively, v is the ion velocity, g the net generation rate, rrii 
the ion mass, V the electric potential, T e the electron 
temperature, and no the electron density at the system’s 
centre. The generation rate g is nil in the inner region I 
but equals gu > 0 in the outer region II. 

Substituting the Boltzmann relation into the other 
equations reduces the system to three equations in three 
unknowns (ji;, v, and V). Doing so, and applying the V 
operator’s definition in cylindrical coordinates, 
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using the fact that, by cylindrical symmetry, v = vr and 
W = (dE/dr)f. 

Supplementing these equations with boundary condi¬ 
tions determines specific solutions. Cylindrical symme¬ 
try dictates that all vector observables, including mv and 
VE, vanish at r = 0. V is zero too at r = 0, for other¬ 
wise the definition of no would be inconsistent with eq. 
©■ The outermost boundary condition comes from the 
demand that solutions be steady states: when the sys¬ 
tem is at a steady state, ions and electrons must leave 
through the outer wall at equal rates, so at r = 2a 


( eV \ pPT7 

„,„ = „ 0 e s p^_jv'^ (8) 

where 2a is the radius of the entire system. At the in¬ 
terface at r = a we require continuity of V (a feature of 
models employing the electron Boltzmann relation) and 
the ion flux niV. 

We nondimensionalize the equations with several char¬ 
acteristic quantities: a generation intensity, go', the Bohm 
speed, c s = ( ksTe/mi ) 1 / 2 ; an ionization length, L; = 
noc s /go; and an electric potential, Vt = fcsT e /e. We 
may then rewrite the above equations in normalized form 
with the dimensionless quantities U = v/c s , G = g/go, 
R = r/Li , A = a/Li, <f> = —V/Vr, and N = ni/no- 


d 2 <F 

dR 2 


d(NU) , NU 
d R + R 

= G 

(9) 

d/7 

U 

d<F 

(10) 

U-^ + G— = 


dR 

N 

dR 

ld$ N - 
+ RdR ~ 

exp(-<F) 

A 2 

(ID 


where A = Xn/Li, Xd being the plasma’s (electron) De¬ 
bye length. The boundary conditions are then 


d$ 

NU = — = $ = 0 

an 


( 12 ) 


at R = 0, continuity of $ and NU at R = A, and at the 
wall at R = 2 A 


NU =^^eM-*) (13) 

The continuity of NU at R = A does not imply con¬ 
tinuity of N or U individually. A priori, N and U may 
have jump discontinuities at R = A. It is similarly im¬ 
portant to note that in the quasi-neutral approximation, 
which takes A —>■ 0 and Hi = n e , eq. (fill) becomes super¬ 
fluous and the wall sheath, being infinitely thin, coincides 
with the wall itself; the R = 2A boundary condition is 
then given by the singularity at which d/7/d R —>■ +oo 
and U —>• 1, the (normalized) Bohm velocity. 


III. SOLUTION FOR REGION I 


The inner region’s solution is trivial because there G = 
0, reducing the continuity and momentum equations to 


and 


d(NU) NU 
d R + IT 


dU _ dd> 
U dR ~ dR 


(14) 

(15) 


respectively. Given eq. m at R = 0, the continuity 
equation’s solution is NU = 0 everywhere in region I, 
and invoking the physical condition that N be every¬ 
where finite implies U = 0 throughout region I. With 
U constantly zero, the momentum equation means $ is 
constant in region I so $ = 0 everywhere in region 
I, not only at R = 0. From the Boltzmann relation we 
then have N = 1 through all of region I. Every variable 
is therefore constant in region I: TV = 1 and U = $ = 0. 


IV. QUASI-NEUTRAL SOLUTION FOR REGION II 


The outer region II, where generation takes place, is 
rather less tractable. There appears to be no closed-form 
solution for eqs. © (EH) here. We may simplify things 
by assuming quasi-neutrality, where we take Hi = n e and 
replace eq. E3 with N = exp(—4)). This leaves two 
equations to solve for N and U. Eq. m, the continuity 
equation, is unchanged but the ion momentum equation 
becomes 


u^L + g Z = _1^L 

dR N N dR 


(16) 


As in region I the continuity equation has a simple 
general solution for NU. With NU continuous at the 
interface of the two regions, region I’s solution implies 
NU = 0 at R = A and the specific solution 


NU = 


G u R 2 - A 2 
~2 R 


(17) 


in region II. The natural choice of go is gu, making Gu = 
1 and 


U 2 R 

Substituting into eq. m, 

(l \dU 2R 2 U 2 + R 2 + A 2 
\U~ JdR~ R(R 2 - A 2 ) 


(18) 


(19) 


solving which gives U as a function of A and R. Again 
there is no closed-form solution, but we can make some 
inferences about the solution’s behaviour. Rearranging, 


dU U 2R 2 U 2 + r 2 + a 2 
dR ~ 1 - U 2 R(R 2 - A 2 ) 
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which is evidently positive where R > A and 0 < U < 1. 
At R = A itself, dU/dR = 1, as shown by rearranging 
eq. m f° r substituting into eq. ([20]) ’s right-hand side, 
and observing that N = exp(O) = 1 at R = A, because 
<f> = 0 in region I and is continuous at R = A. 

As such U strictly monotonically increases with R until 
U = 1, at which point dU/dR —>• +00 and there is a 
singularity. We thus reproduce the usual quasi-neutral 
behaviour: the ion speed tends to the Bohm speech as 
the electric field tends to infinity. 

For want of an analytic solution we solve eq. m nu¬ 
merically. To accelerate the numerical computation we 
exploit the fact that given U = Uq at some position, eq. 
m may be approximated as 

(1 TT \ dU 2R 2 U 2 + R 2 + A 2 
\U Jd R R(R 2 — A 2 ) 1 1 

where U remains a function of R, but Uq is constant. 
Applying the condition U = Uq at R = Rq, this (now 
separable) differential equation has the solution 


U = 


\ 


~ W \ exp ( _C/ o) 


^42 _ R'2\ 2 ( 1 + !7 o) 


A 2 -R 2 


( 22 ) 

where W is the Lambert W function’s principal branch. 
This solution is valid for any interval [Rq , R] over which 
U ~Uq. We may then obtain the quasi-neutral solution 
over all of region II by chopping region II into subregions, 
all so narrow that U is virtually constant within each, and 
applying eq. (E^l) to each subregion in turn to compute 


U. 


We simultaneously compute N from eq. m, and $ as 
In (1 /TV). Given the boundary condition that U = 1 at 
R = 2A, the choice of A determines N and hence $ at 
the interface. Since $ = 0 there, we must then choose 
A to make IV = 1 at R = A. By numerical trial and 
improvement we find this critical A value is 0.62542, so 
the total radius of the quasi-neutral system is 1.251Ti, 
and N at its outermost edge is 0.468. Figure [2] presents 
this canonical solution for U and N. 

The quasi-neutral solution does not hold near the sys¬ 
tem’s outer wall, where the plasma approximation fails - 
in reality a sheath must form there, so rij ceases to be ap¬ 
proximately n e . However, in the A — » 0 limit the sheath 
is vanishingly thin and the quasi-neutral solution gives 
excellent results almost everywhere. The quasi-neutral 
solution therefore provides a rough check on numerical 
solutions with small A. 


V. NUMERICAL SOLUTION FOR REGION II WITH 
FINITE DEBYE LENGTH 

Because A is nonzero in practice, realistic solutions of 
eqs. ® ED differ from the quasi-neutral solutions, and 
the equations must be re-solved numerically. To accom¬ 
plish this we apply the shooting method. In each itera¬ 


tion of the shooting method we solve the three differential 
equations by integrating outwards from near the interface 
with the midpoint method. 

The normalized system has the four dimensionless pa¬ 
rameters mi/m e , On, A, and A; once any three are set 
the fourth is then constrained to a particular value by the 
outer wall boundary condition. We fix mi/m e = 736744 
(corresponding to Hg(/ ions, the ion species considered 
in our previous papers), Gu = 1, and A = 0.003, then 
find A by trial and improvement, converging to an A 
for which eq. (1131) is approximately satisfied. Figure [3] 
presents this solution of the system. We find that the 
region II plasma has a sheath adjacent to the outer wall, 
but has minimal net space charge closer to the plasma’s 
central axis. Near the interface (R ss A) the numerical 
solution is similar to the quasi-neutral solution, but the 
two solutions gradually deviate farther out. 

A subtle but important feature is a discontinuity in 
N at the interface. In region I, N is exactly unity. 
However, for our current parameters, N discontinueusly 
jumps from 1 to 1.000018 at the beginning of region II. 
The electron density is continuous at the interface be¬ 
cause it is given by the Boltzmann relation, and we can¬ 
not have discontinuities of potential. This argument does 
not apply to the cold positive ions. The ion density must 
change in order to produce a very small space charge 
density (a divergence of the electric field) in the outer 
plasma. We can note that this does not imply a layer of 
charge on the boundary, the electric field is continuous. 
Consequently, all of region II is a positive space charge 
region, not just its sheath, but that space charge is too 
subtle to be visible near the interface unless one looks 
closely (figure 03 inset). 

TV’s discontinuity is not an artifact. Even under quasi¬ 
neutrality, N and $ have nonzero curvature at the in¬ 
terface, which in physical terms implies a nonzero space 
charge density. To estimate that density, we rearrange 
eq. ED for U and substitute into the quasi-neutral ion 
momentum equation. After some algebra that gives 


— IV 


dTV 

dR 


r2 _ A 2 R 2_ A 2 

2 R + 4 R 3 N 


(R 2 + A 2 ) N 


- (R 2 - A 2 ) R 



(23) 


At R = A, dTV/dR = 0 and R 2 — A 2 = 0, so at that R 
all terms in eq. (l23l) are nil. But the terms do not all 
have the same dependence on R — A. Unlike the other 
terms, the very last term is a product of dN/dR and of 
(R 2 — A 2 ) squared, so for 0 < (R — A) -C 1 it is negligible 
relative to the other terms. We therefore drop this term 
to get 


-TV 


dTV _ R 2 -A 2 
d R ~ 2R 


R 4 - A 4 
4 R 3 


(24) 


2 R 
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FIG. 2. Numerical quasi-neutral solutions for U (left panel, solid curve), $ (left panel, dashed curve), and N (right panel) with 
A = 0.6254204. 


which has the specific solution 


N 


4 (Nq + A 2 ) - 3R 2 - ^ + 4A 2 In 


R 

A 


(25) 


where No is N’s value at the interface (R = A). Under 
quasi-neutrality, $ = — In IV, so we substitute the N 
solution into that expression for $ and take $’s Taylor 
series about R = A to second order: 


$ 

d 2 d> 

dR 2 


— In N 0 + 


{R - A) 2 
N 2 



(26) 

(27) 


For all A > 0, v and hence N 0 are greater than 1. It 
follows that N jumps discontinuously from 1 to No where 
region I becomes region II, and the size of the jump is the 
net space charge density at the interface. (Our numerical 
solver uses this fact to compute boundary conditions for 
the region II equations.) When A_d <C Li, N 0 ~ 1 + 2A 2 , 
so the net space charge density at the interface is ss 2A 2 
(figure 0. 

It is straightforward to consider the solution’s spatial 
dependence near the interface rather than its A depen¬ 
dence. Taking Taylor series of eq. (E5l) about R = A, 

N = N 0 - ( R ~ A ' ) + O ({R- A) 3 ) (32) 

No 


<!>’s curvature is not going to be much different for suffi¬ 
ciently small A, so we now let A be finite and reintroduce 
Poisson’s equation: 


d 2 <4> _ N - exp(—<h) 1 d$ _ 2 
dR 2 ~ A 2 R dR ~ N$ 


(28) 


At the interface, N = N 0 and = d$/dR = 0, so 


No-1 ^ 2 
A 2 ~ N 2 

a cubic equation with the solution 


(29) 


where 


iV 0 


1 + v + l/v 
3 


(30) 


v = \j 1 + 3Aa/6”+ 81A 2 + 27A 2 


(31) 


In the limit of A vanishing outright, Nq —> 1. We then 
have N ~ 1 — ( R— A ) 2 , $ ~ ( R — A ) 2 from eq. (EBll . and, 
combining the continuity equation with the fact that U = 
0 and G jumps to 1 at the interface, dU/dR « 1/N ss 1, 
so U ~ (R—A). This matches the behaviour of the curves 
in figurc[2]just beyond the interface: N falls parabolically, 
<f> grows parabolically, and U rises nearly linearly. 

Just as the quasi-neutral case casts light on systems 
with 0 < A <C 1, our numerical solutions of the latter 
support our quasi-neutral result. Computing A for as¬ 
sorted small values of A (with G and rrii/m e fixed at 1 
and 736744 respectively), it’s evident that as A shrinks, 
the resulting A values converge on 0.62542 (figure [5]), 
the canonical A value found earlier by solving the quasi¬ 
neutral equations. For A <C 1, A is greater than the 
canonical quasi-neutral A by about 4A 4 / 5 . 

Figure [5] displays the same information in a different 
manner. The quantity 1/A illustrates that the ionization 
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R = r/L, 
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FIG. 3. Numerical region II solution for A = 0.003 (thin curves) compared to the canonical quasi-neutral solution (thick dotted 
curves). We augment the default choice of axis scale (R, on the lower axes) with the alternative choice (on the upper axes) of 
r normalized by a’s value when A = 0.003. In plot (b), the thin solid curve represents the ion density while the thin dashed 
curve represents the electron density exp(—4>). 


length is comparable to the plasma dimension, but de¬ 
pends on Ad’s magnitude in comparison with the latter. 


VI. CONCLUSION 

An interesting result of the present investigation is 
that a uniform field-free plasma forms in the inner re¬ 


gion. This is an extremely rare occurrence in plasma 
physics, the other example known to us being that of a 
thermally produced caesium plasma in a cavity^. It is 
clear that further calculations could readily be carried 
out. An example would be to use a Tonks-Langmuir 
model^ in which the positive ions have a distribution of 
velocities in the radial direction. Collisions with neutral 
atoms could also be included in the analysis. A process 
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FIG. 4. Detail of the ion density discontinuity at the interface 
between regions I and II in the A = 0.003 case with HgJ ions 
and G = 1. 



not included in the model is that electrons produced by 
irradiation may then contribute to further ionization in 
collisions with neutral atoms. 

We have also confined ourselves in this paper to the 
system’s steady-state solutions, but one could carry out 
a transient analysis which accounted for the initial flow of 



The dashed line is 1/A in the A —> 0 limit. 

electrons into region I as the plasma establishes equilib¬ 
rium. This electron flow would produce an electric field 
drawing the ions in after the electrons, a phenomenon 
akin to the expansion of a plasma into a vacuum^. 
Modelling transient effects could also quantify how the 
plasma’s behaviour adjusts the fundamental, observable 
plasma parameters (no, T e , and gu) to move A to its 
steady-state value, even if the experimenter’s initial con¬ 
ditions give A a different starting value. 

Calculations could be carried out for a less abrupt 
change between regions I and II as briefly discussed in 
our earlier paper—, and for different ratios of the two re¬ 
gions’ radii. We made our choice of positive ion, i.e. HgJ, 
for historical reasons given in the same paper—. The ratio 
of this ion mass to the electron mass enters the theory 
through the boundary condition at the wall. In effect it 
determines only the width of the sheath and the voltage 
across it. It is of interest to note that the ion velocity on 
entering the sheath is about equal to the Bohm velocity, 
whereas in our earlier model, which had ionization in the 
inner region, considerably higher velocities were foundA. 
An even more striking difference between that model and 
the current model is the behaviour of their solutions near 
the interface. Our earlier model produced a double layer 
centred on the interface, while our current model has a 
jump discontinuity in density. 

Our previous paper gave an indication of the equip¬ 
ment needed for experimental work in this field-. The 
model system we present here should be relatively easy 
to realize in the laboratory, especially compared to our 
original model system, and could be a way to make uni¬ 
form, stationary plasmas on command. To accomplish 
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this, annular UV illumination of a cylindrical container 
of mercury gas should suffice. We are hopeful that both 
experimental work and further theoretical work shall be 
carried out to make further progress in this vein. 
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